Asynchronous entanglement from coherently coupled nonlinear cavities. 
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The output fields of a pair of coherently coupled nonlinear optical cavities are found to exhibit strong opti- 
cal entanglement. For sufficiently strong coupling the quantum correlations become asynchronous providing a 
resource for quantum information protocols such as all-optical quantum memories. A straightforward experi- 
mental implementation applicable to whispering gallery mode resonators such as microtoroids is proposed. 



Last century, interest in entanglement was motivated pre- 
dominately by its fundamental significance as "the character- 
istic trait of quantum mechanics"! 1 ]. In recent years, however, 
a paradigm shift has occurred due to the realization that entan- 
glement can facilitate powerful measurement, computation, 
and communication tasks[2|. Optical entanglement is of par- 
ticular significance, both enabling applications such as quan- 
tum cryptography, quantum information networks, and quan- 
tum metrology[2|; and facilitating the most stringent tests of 
quantum mechanics performed to date 0. 

Here we investigate, for the first time, optical entanglement 
generated from a pair of coherently coupled nonlinear opti- 
cal cavities. Coherently coupled optical cavities are of fun- 
damental interest, enabling for example an all-optical analog 
to electromagnetically induced transparency (EIT)(4l with the 
capacity to both slow and stop lightQ. We show that the 
introduction of nonlinearity enables arbitrarily strong entan- 
glement to be generated. Interestingly, however, for strong 
coupling, the entanglement is quenched, and a new form of 
asynchronous entanglement exhibiting time-delayed quantum 
correlations becomes apparent. This entanglement is analo- 
gous to that generated in A-type atomic ensemble quantum 
memories! 6], which have already proved to be an enabling 
technology in quantum information applications such as quan- 
tum repeaters and on-demand single photon sources 171. 

In contrast to atomic ensembles, quantum memory capabil- 
ities are not typically available in nonlinear optical systems. 
Indeed, they are only present in our system due to the EIT- 
like nature of the coupled cavities; and to our knowledge have 
not been predicted for any other nonlinear optics based entan- 
glement source. The asynchronous entanglement predicted 
here has many applications in quantum information science. 
For example, an all-optical quantum memory could be im- 
plemented utilizing asynchronous entanglement as the non- 
classical resource for quantum teleportation[8|. 

This Letter models a pair of identical below threshold co- 
herently coupled optical parametric oscillators (OPOs) each 
consisting of a ^ k+1 ) general nonlinear medium enclosed in 
an optical cavity as shown in Fig.[TJa). Such a configuration 
is experimentally relevant, most notably to ultrahigh quality 
whispering gallery mode (WGM) resonators where counter 
propagating modes are coupled by scattering centers as 
shown in Fig. [TJb). Silica microtoroidal resonators [9], and 
polished crystal microdisks fT3lL for example respectively ex- 




FIG. 1: (color online). Schematics of coupled cavity systems, a) two 
optical cavities coupled via a partially reflective mirror, each con- 
taining a (/c + l) th order nonlinear medium, b) an equivalent WGM 
resonator with two counter propagating optical modes coupled by 
optical scattering and intrinsic nonlinearity. 

hibit high x^ and x^ nonlinearity, as well as strong scat- 
tering. These resonators are capable of providing strong op- 
tical confinement in a scalable, robust, microfabricated ar- 
chitecture with high coupling efficiency to optical fiber|9ll. 
They are therefore an ideal candidate to produce, in-fiber, 
the strong entanglement predicted here. Such entanglement 
would be an important resource for scalable quantum infor- 
mation networks, as well as long-distance fundamental tests 
of quantum mechanics such as the EPR paradox ifTUl and 
Bell testsQ. Furthermore, recent research has demonstrated 
an exquisite capability to engineer microresonator mechani- 
cal properties ifTTIl . providing the potential to suppress noise 
due to guided wave Brillouin and Raman scattering which 
presents a strong constraint on existing in-fiber entanglement 
sources fl2l . 

Analysis: In our model, within each coupled OPO the non- 
linear process converts k pump photons at frequency uo p into a 
pair of signal and idler photons at frequencies u\ = ft— A' and 
002 = ^+A'; where to maintain energy conservation Q = |cOp. 
Since the pump fields in such systems are typically bright we 
treat them classically here, with a x and a y denoting their co- 
herent amplitudes, and the subscripts x and y used throughout 
to distinguish the two cavity modes. Assuming the nonlinear 
interaction strength F is identical for modes x and y, and the 
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coherent coupling rate g is independent of detuning A', as is 
the experimentally relevant case, and taking the rotating wave 
approximation yields the system Hamiltonian 



g(a xn a yn +a xn a yn ) + ^ ^ uj rnn d^ mn a ri 
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which combines the Hamiltonians describing QPQ 11141 [151 
and coherent coupling lfT6l used in previous work; where m G 
{x, y}, a is the annihilation operator, and the subscripts 1 and 
2 denote the signal and idler fields respectively. Assuming 
both cavity modes x and y have identical input coupling and 
loss rates, denoted by 7i n and respectively, as is typical 
for experimental systems, and applying the quantum Langevin 
equation ifTTl to Eq. ([T]) we obtain four equations of motion for 
the signal and idler fields 

a xl = Ta k x a x2 - (j+iA)a xl - ig'a yl - y^Ti^Wi + >ffiicn iXl 
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a X 2 = Ta*a} xl - (~/-iA)a x2 - ig'a y 2 ~ V^Th^in 
Oyi = Tay]a y2 - (7+zA / )a 2/ i - ig f a xl - v27i^a in , y i + yfeyiai iy i 
a y2 = rajaji - (j-iA)a y2 - ig'd x2 - v27i^a in , y2 + \finia,i iV 2 

where a- in and di are the fields entering through the input 
coupler and loss mechanisms respectively, and 7 = 7i n + 7/ 
is the total decay rate of each cavity. Assuming a x and a y 
are independent of time these equations can be solved in the 
frequency domain by taking the Fourier transform. 

The dimensionless quantities R x e lk< ^ x = r l^l e ^<j> x ^ 

are used; where 77 is the cavity escape efficiency and R is a di- 
mensionless pump strength with the physical significance that 
R = 1 corresponds to the OPO threshold. The solution can 
then be expressed compactly as \MA = y/rjA{ n + y/l — rjAi, 
where 
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and throughout vectors A q are of the form A q = 

[ a q,xl, a q yX 2i a <?,2/l' a q,y2\ * 

One can then use the input/output relations IfTTl A out = 
v /27 in A+A in to solve for the output fields, 



A out = [h + 2r ? M- 1 ]A in + 2^r ] (l-r ] )M- 1 A h (3) 

where 1 4 is the 4x4 identity matrix. General output field 
quadratures for each of the four system modes are then readily 
obtained as X 9 (A) = a( A)e~ ie +a f ( A)e iG . 



Gaussian continuous variable entanglement between or- 
thogonal quadratures pairs of two light beams can be 
completely characterized by the correlation matrix of the 
system lfT4l ; which can be constructed in our case from Eq. ([3} 
for any two pairs of output fields. Here we follow the proce- 
dure of Duan et. al [18 ] for bipartite continuous variable sys- 
tems, who showed that the correlation matrix of any Gaussian 
state can be transformed reversibly into a standard form, from 
which the degree of inseparability X can be defined such that 
X < 1 is a necessary and sufficient condition for entanglement, 
with maximal entanglement characterized by X = 0. Here we 
follow the method of Grosse| 15 ], to numerically calculate the 
degree of inseparability from the output field correlation ma- 
trix of the coupled cavity system. 

Results: In general, entanglement can be observed between 
the fields a x i, a x2 , a y i, and a y2 through direct individual 
measurements on each. However, here we consider the more 
experimentally realistic situation of heterodyne detection per- 
formed separately on the outputs of each of the subsystems x 
and y with a local oscillator at frequency 

In this case for each subsystem both signal and idler fields 
a m i and a m2 beat with the local oscillator and, combined, 
produce a signal directly proportional to the joint-quadrature 
operator 



x e m AA)= 1 
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where m G y}, is the local oscillator phase, and t is 
the time at which the measurement occurs. As a result of 
the frequency difference between the signal and idler fields 
we see that a time delay has the effect of rotating the ob- 
served signal and idler quadratures in opposite directions. En- 
tanglement of the orthogonal quadrature pair X^ x (A) and 

X x x ^^ 2 (A) of field x with the equivalent quadratures of 
field y is possible due to the Heisenberg uncertainty relation 
A 2 X^ n A 2 X^~ 7V ^ 2 > 1 which originates from the commuta- 



tion relation [Xf^ , X 



0+7r/2i _ 



2i. We analyze entanglement of 



this kind henceforth, and restrict our analysis to the optimal 
situation where x =0 y =0. 

In all nonlinear optics based entanglement sources to date, 
due to the near instantaneous response of the nonlinearity the 
quantum correlations between output fields are observed to 
be synchronous. Hence, here we begin by considering the 
case where the output fields are measured at the same time, 
i.e. t x = t y . Fig. [2ja,c) illustrates the degree of insepara- 
bility calculated for this case as a function of coupling rate 
and detuning from cavity resonance, where we have chosen 
equal intensity in-phase pump fields and parameters typically 
achieved in experiments. 

Strong entanglement is observed. When g <C 1 the best 
entanglement is near resonance as one might perhaps expect. 
However for large coupling one observes both that the en- 
tanglement is quenched and that the maximum entanglement 
shifts to nonzero detuning. The shift can be understood from 
the well known mode splitting that occurs in coupled cavity 
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FIG. 2: (color online). Inseparability as a function of detuning A 
and coupling rate g. a) and c) synchronous entanglement with Acp = 
0; b) and d) asynchronous entanglement with Acj> = ir/k. Model 
parameters: R x — R y — 0.9 and r\ — 0.99. 




Time delay tx[A/ti] 

FIG. 3: (color online). Inseparability as a function of measurement 
time delay optimized over detuning, a) Acj) = 0, b) Acj> = 7v/2k, c) 
A<j) — n/k. Model parameters: R x — R y — {).9, r/ = 0.99, g = l. 

systems (9). As a result fields produced at frequencies de- 
tuned by g from the unperturbed resonance frequencies are 
on resonance within the cavity and the nonlinear response is 
enhanced. The entanglement quenching, however, is surpris- 
ing since the coupling is a reversible process. To investigate 
this further we consider the asynchronous case where the mea- 
surement on field x is delayed by a time r = t x —t y . 

Fig. [3] shows the evolution of entanglement with measure- 
ment time delay for three different relative pump phases Acj), 
where we have optimized the inseparability over detuning A 
as is typically performed in experiments using spectral analy- 
sis. In Fig.[3ja) the pumps are in phase and we see the normal 
behavior for a nonlinear optics based entanglement source, 
with maximal entanglement in the synchronous case where 
r = 0; and as seen before in Fig.[2ja,c) a strong dependence 
on g, with entanglement optimized at g « 0.5. In Fig.JSJb) a 
pump phase delay of tt/ 2 k is introduced and we see an over- 
all reduction in the level of entanglement. However, now as 
the coupling rate increases and the synchronous entanglement 
is quenched, a transition occurs to entanglement in an asyn- 
chronous regime with a time delay between detection events 
of r = 7r/2A. In contrast, this asynchronous entanglement is 
optimized at large g. In Fig.J^c) where Acj) = n/k the syn- 
chronous entanglement is no longer apparent for any g, with 
asynchronous entanglement present for all but very small g. 
The dependence of the degree of asynchronous inseparabil- 
ity on coupling rate and detuning for this case are shown in 



Fig. [2jb,d). It is clear that now the maximum entanglement 
always occurs at A = g, can be shifted far from the resonance 
frequency of the cavity, and is insensitive to g for g ^> 1. This 
resonance frequency shift is technically significant, allowing 
the optimal entanglement to be moved away from low fre- 
quency noise sources such as the laser relaxation oscillation 
which typically limit both entanglement strength and purity. 

Fig. [4] shows the dependence of the degree of inseparabil- 
ity on pump strength and cavity escape efficiency at optimal 
A. For both synchronous and asynchronous cases, we choose 
the coupling rate to provide maximal entanglement, and see 
that complete entanglement is achieved as the pump strength 
and escape efficiency approach unity, although with different 
functional dependance. This functional dependence will be 
studied further in future work. 

A standard method of generating optical entanglement is to 
combine a pair of squeezed fields on a 50/50 beam splitter |T9l , 
with the relative phase being critical to the strength of the en- 
tanglement. In an analogous manner the coherent coupling in 
our scheme couples the non-classical fields generated in the 
two cavity modes. It is therefore to be expected that the rel- 
ative pump phase Acj) will critically affect the entanglement 
achievable in this system. Fig. [5] shows the minimum degree 
of inseparability between the output fields optimized over time 
delay and detuning as a function of relative pump phase and 
pump strength for different orders of nonlinearity. A cou- 
pling rate of g = 1 is chosen which, as shown in Fig. |3j is 
an interesting case where strong entanglement can be simul- 
taneously generated both synchronously and asynchronously. 
This represents a regime of asynchronous 4-partite entangle- 
ment, where the time-delayed and the non-time-delayed out- 
puts of x are both entangled to both time-delayed and non- 
time-delayed outputs of y. We see that 4-partite entanglement 
is apparent throughout the majority of the parameter space. 
This multipartite entanglement will be investigated further in 
future work. 

Returning to bi-partite entanglement, it is clear from 
Fig. [5] that the entanglement increases with increasing pump 
strength. We also see a strong dependence on A <j> as expected, 
however for high pump strength substantial entanglement can 
be generated for all Acj) in distinct contrast to conventional en- 
tanglement sources where typically entanglement generation 
fails over some range of Acj). 

Finally we turn our attention to a robust and straightforward 




FIG. 4: (color online). Inseparability optimized over detuning as a 
function of pump strength R — R x — R y and cavity escape effi- 
ciency r\. a) synchronous entanglement with g = 0.5 and Acj) — 0, 
b) asynchronous entanglement with g ^> 1 and Acj> = ir/k. 
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FIG. 5: (color online). Dependence of inseparability on pump strength R = R x = R y (radial axis) and relative phase A</> (azimuthal axis), 
optimized over A and r. Order of nonlinearity: a) %^ (k = 1), b) (k = 2), c) x^ (k = 3). Dark lines circumscribe sectors of 4-partite 
entanglement; grey bands circumscribe sectors dominated by asynchronous entanglement. Arrow: path followed for single sided pumping. 
Model parameters: 77 = 0.99 and g— 1. The plot sectors are effectively compressed by ir/k as the order of the nonlinearity increases. 



experimental implementation. Consider the case where only 
cavity mode x is directly pumped with mode y pumped indi- 
rectly through coherent coupling from mode x. In this case 
one finds that a y = iga x . Hence the relative phase between 
the pumps is A(j) = 37r/4, so that when g = 1 single sided 
pumping tracks the downwards facing arrow in Fig. [5] We 
see that in the special case of x^ nonlinearity relevant to sil- 
ica microtoroidal resonators [ 9 ] optimum entanglement can be 
achieved naturally in this single sided pumping scenario. For 
this case the simple analytical solution 




16r]R[A - R(l - 77)] 

[(1 + RY + (A - 1) 2 ][(1 - RY + (A + l) 2 ] ' 



for the inseparability can be obtained. It is clear that as 
{i2,77,A}— >0. Hence for sufficiently strong pump- 
ing and high efficiency, strong entanglement can be generated 
with only a single external pump field. 

In conclusion, we have shown that strong entanglement can 
be generated between the output fields of a pair of coherently 
coupled nonlinear optical cavities. For sufficiently strong cou- 
pling, the quantum correlations become asynchronous, and 
can therefore be used in conjunction with quantum telepor- 
tation to achieve an all-optical quantum memory. A straight- 
forward experimental implementation applicable to ultrahigh 
quality WGM microresonators such as silica microtoroids is 
proposed requiring only one external pump field. 
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